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ABSTRACT 

The statistical expectation values of the temperature fluctuations of the cosmic microwave background (CMB) 
are assumed to be preserved under rotations of the sky. We use the bipolar power s pectrum (BiPS) proposed in 
Hajian & Souradeep to investigate the statistical isotropy (SI) of the CMB anisotropy maps recently measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP). The method can isolate and probe specific regions of choice 
in multipole space using appropriate window functions. The BiPS is estimated for full sky CMB anisotropy maps 
based on the first year WMAP data using a range of window functions. The BiPS spectra computed for both full 
sky maps for all our window functions are consistent with zero, roughly within 2 a. The null BiPS results may 
be interpreted as an absence of strong violation of statistical isotropy in the first-year WMAP data on angular 
scales larger than that corresponding to / ^ 60. However, pending a careful direct comparison, our results do not 
necessarily conflict with the specific SI related anomalies reported using other statistical tests. 

Subject headings: cosmic microwave background - cosmology: observations 



1. INTRODUCTION 

The cosmic microwave background (CMB) anisotropy is 
a very powerful observational probe of cosmology. The re- 
cent Wilkinson Microwave Anisotropy Probe (WMAP) data 
have provided a precise determination of standard cos- 
mological parameters reaffirming an emerging concordance 
model of cosmology, structure formation (Bennett et al. 2003; 
Spergel et al. 2003) and, even inflation and the early universe. 
This remarkable success story is almost entirely based on the 
measurements of the angular power spectrum, C/, of the CMB 
anisotropy. 

In standard cosmology, the CMB anisotropy Ar(n) is ex- 
pected to be statistically isotropic, i.e., the statistical expecta- 
tion values of AT(fi) are preserved under rotations of the sky. 
The issue of statistical isotropy (SI) applies to both Gaussian 
and non-Gaussian CMB sky maps. Here we limit our study to 
the Gaussian CMB anisotropy field, where the two-point corre- 
lation function contains all the statistical information encoded 
in the field. The importance of SI is based on a very well known 
fact: if the two point correlation is rotationally invariant, one 
can equivalently completely specify CMB anisotropy in terms 
of the Legendre transform, C/, the widely used angular power 
spectrum of CMB. When the SI holds, the information of the 
full CMB data set can be compressed into a multipole spec- 
trum. 

After the recent release of the first year of WMAP data, the 
SI of the CMB anisotropy has attracted considerable atten- 
tion. Tantalizing evidence of SI breakdown (albeit, in very 
different guises) has mounted in the WMAP first year sky 
maps, using a variety of different statistics. It was pointed out 
that the suppression of power in the quadrupole and octopole 
are aligned (Tegmark et al. 2004). Further "multipole-vector" 
directions associated with these multipoles (and some other 
low multipoles as well) appear to be anomalously correlated 



(Copi et al. 2004; Schwarz et al. 2004). There are indications 
of asymmetry in the power spectrum at low multipoles in op- 
posite hemispheres (Eriksen et al. 2004a; Hansen et al. 2004; 
Naselsky et al. 2004). Possibly related, are the results of tests 
of Gaussianity that show asymmetry in the amplitude of the 
measured genus amplitude (at about 2 to 3(t significance) be- 
tween the north and south galactic hemispheres (Park 2004; 
Eriksen et al. 2004b). Analysis of the distribution of extrema 
in WMAP sky maps has indicated non-gaussianity, and to some 
extent, violation of SI (Larson & Wandelt 2004). However, 
what is missing is a common, well defined, mathematical lan- 
guage to quantify SI (as distinct from non Gaussianity) and the 
ability to ascribe statistical significance to the anomalies unam- 
biguously. 

These results provide considerably increased motivation 
for this Letter in which we test the WMAP data for SI, in 
detail, using the bipolar power spectrum (BiPS) of CMB 
anisotropy, a powerful and easy to understand statistics, 
which was recently introduced to the field of CMB analysis 
(Hajian & Souradeep 2003b; Souradeep & Hajian 2003). If SI 
is violated, the C/ do not contain all the statistical information 
of a Gaussian sky map. It is described by the whole set of cross- 
correlations of the coefficients in the harmonic space (a/,„fl;/„,'). 
The cosmic variance of the elements will prevent us from esti- 
mating the individual elements from a single map. Hence, it is 
very important to define a statistic that condenses the informa- 
tion of that matrix into a single spectrum that reduces the cos- 
mic variance. BiPS does exactly that. Furthermore it is easy to 
interpret and straightforward to implement for very rapid com- 
putation. 

2. BIPOLAR POWER SPECTRUM ANALYSIS 

In spherical harmonic space representation where AT(n) = 
J^im^imYimin), the condition for SI translates to a diagonal 
{ai,„ap^,) = Ci6ii'S,„m' where C/, is a complete description of 
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(Gaussian) CMB anisotropy. When {ai„ap^,) is not diagonal, 
the off-diagonal elements contain information beyond Q . 

The BiPS is a combination of these off-diagonal elements 
in (aimap^,), containing extra non-SI information. The BiPS 
together with C/ give a complete, orientation independent de- 
scription of a Gaussian CMB anisotropy map. 

The BiPS of CMB is defined by (Hajian & Souradeep 2003b; 
Souradeep & Hajian 2003; Hajian & Souradeep 2004) 



ll'M 



(1) 



where Ajjf are the coefficients of the bipolar spherical harmonic 
(BiPoSH) expansion of the correlation function, which is the 
most general expansion of C(q, q'): 



Ciq, q') = Y, A'lffiYiiq) ® Yv{q')}iM ■ 
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BiPoSH functions are the natural complete basis on S^ x S^ 
(Varshalovich et al. 1988). Here A^iff are a complete set of lin- 
ear combinations of off-diagonal elements of {aimO*,^,): 



A« = ^(fl,„flrw)(-i)"''e:; 
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where €.^i^i,^, are Clebsch-Gordan coefficients. Every spe- 
cific measure of off-diagonal elements of the covariance ma- 
trix {aimtt*,^,) suggested in the literature (Durreretal. 1998; 
Prunet et al. 2004) can be expressed in terms of Ajff . For 
an SI correlation function {aimOp^,) = Cidu'dmm', implying 
A« = {-iyCii2l+iy^Hii, SmSMO and Ke = n^m- An equiv- 
alent real space description of the BiPS in terms of ap- 
propriately weighted averages of the correlation C(q,q') is 
given in previous publications (Hajian & Souradeep 2003b; 
Souradeep & Hajian 2003). 

We compute the BiPS using fast methods of spherical har- 
monic transform of the map ' and apply filters in multipole 
space using the appropriate positive definite window function 
Wi. We define an unbiased estimator for the BiPoSH coeffi- 
cients and then estimate k^, 



4f = VM^^fltaflr„. c; 
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have restricted our analysis to Z < 60 where the errors are dom- 
inated by the cosmic variance. 

We carry out measurement of the BiPS, «£, on two full sky 
CMB anisotropy maps - (A) the Internal Linear Combination 
map (denoted as "ILC" in the figures; (Bennett et al. 2003)), 
and (B) a foreground cleaned map (denoted as "Tegmark"; 
(Tegmark et al. 2004)). The angular power spectra of these 
maps, shown in Figure 1, are consistent with the best fit the- 
oretical power spectrum from the WMAP analysis. ^ 

We use simple filter functions in / space to isolate different 
ranges of angular scales. A low pass, Gaussian filter 



Wi'^ih) = exp[-(Z+ l/2)7(;,+ 1/2)2], 



(5) 



3£ = [Kjl^j, is the bias that arises from the SI part of 
the correlation function. Although BiPS is quartic in a;,,,, it is 
designed to detect SI violation and not non-Gaussianity. 

Consequently, for SI correlation, the measured k£ will be 
consistent with zero within the cosmic variance. Cosmic er- 
ror, a^^, and bias, Sf, are given by equations (14) and (17) 
of (Hajian & Souradeep 2003b) respectively, with Q -^ WiCi, 
where Ci = Cf + Cf is the total SI angular power spectrum of the 
"true" signal and noise. Multipole space windows that weigh 
down the contribution from the rest of the SI region of the mul- 
tipole space will enhance the signal relative to cosmic error, a^^ . 

It is important to note that bias can never be subtracted ex- 
actly for a non-SI map. What is important is whether the mea- 
sured Ki differs from zero at a statistically significant level. 

It is not inconceivable that for strong SI violation, *B£ over- 
corrects for the bias leading to negative values of kg. The noise 
covariance is a possible source of SI violation. However, we 

1 We use the ANAFAST routine of the HEALPix package (Gorski.Hivon & Wandelt 
^ Based on an LCDM model with a scale-dependent (running) spectral index. 



which cuts off power on small angular scales (< l/ls), and a 
band pass filter, 

W,'(h,l,) = (2{l-M(l+l/2)/(h+ 1/2)]})W,«(« (6) 
which retains power within a range of multipoles set by It and 
Is- The window functions used in our work are plotted in Figure 
1. 

We use the WMAP best-fit (WMAPbf) Q to generate 1000 
simulations of the SI CMB maps. The «/„, are generated up 
to an I of 1024 (corresponding to maps at HEALPix resolution 
Nside = 512). These are then multiplied by the window func- 
tions Wf'(ls) and Wi^(lt,ls) and the BiPS for each realization is 
computed. We use Cf to analytically compute the bias and the 
cosmic variance for Kg. The average and standard deviation of 
Kf of the SI maps are also estimations of the bias and the cos- 
mic variance. We verified that the theoretical cosmic variance 
and bias match the numerical estimations of standard deviation 
and average Kg of the 1000 realizations of the SI maps. 

The BiPS of the full-sky CMB maps (A) and (B) based on the 
WMAP first-year data is computed in the same manner. We use 
the WMAPbf Cf to compute the bias and cosmic variance ana- 
lytically. Using the analytical bias and cosmic variance allows 
us to rapidly compute BiPS with 1 a error bars for different the- 
oretical Cf . Figure 2 shows the measured values of Kg for maps 
(A) and (B) for two of the window functions. We compute the 
BiPS for all the window functions shown in Figure 1 . 

The BiPS measured from the 1000 simulated SI realizations 
of WMAPbf Ci is also used to estimate the probability distri- 
bution functions (PDFs), defined as p{ki) = /- dnipini) for 

K; > and j'^' dKip(Ki) for «/ < (Fig. 3). We compute the 
individual probabilities of the map being SI for each of the mea- 
sured Kg. 

The probabilities for the W;'^(20,30) window function are 
greater than 0.25 and the minimum probability at ^ 0.05 oc- 
curs at K4 for W;'-'(40). The reason for systematically lower SI 
probabilities for W/'^(20,30) as compared to W/'^(40) is simply 
due to the lower cosmic variance of the former. The contribu- 
tion to the cosmic variance of BiPS is dominated by the low 
spherical harmonic multipoles. Filters that suppress the «/„, at 
low multipoles have a lower cosmic variance. 

It is important to note that this probability is a conditional 
probability of the measured kg being SI given the theoretical 
spectrum Cf used to estimate the bias. A final probability 
emerges as the Bayesian chain product with the probability of 
the theoretical Cf used given data. Hence, small differences 
in these conditional probabilities for the two maps are perhaps 
not significant. The important role played by the choice of the 
theoretical model for the BiPS measurement is shown for a W/ 

1999) publicly available at http://www.eso.org/science/healpix/. 
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that retains power in the lowest multipoles, I = 2 and 3. Assum- 
ing WMAPbf Cj , there are hints of non-SI detections in the 
low I (Fig. 4, top). We also compute the BiPS using a Cj for 
a model that accounts for the suppressed quadrupole and oc- 
topole in the WMAP data (Shafieloo & Souradeep 2004). The 
mild detections of a non zero BiPS vanish for this case (Fig. 4, 
bottom). This suggests that the excess power in the WMAPbf 
Cj with respect to the measured Q from WMAP at the lowest 
multipoles tends to indicate mild deviations from SI. 

The SI of the CMB anisotropy has been under scrutiny after 
the release of the first year of WMAP data. Using a BiPS anal- 
ysis we find no strong evidence for SI violation in the WMAP 
CMB anisotropy maps. We have verified that our null results 
are consistent with measurements on simulated SI maps. 

It is also possible to construct a model-independent "frequen- 
tist" estimator of BiPS that uses the C; of the map itself to com- 
pute the bias. Measurements on our simulated sky maps show 
that cosmic variance for this estimator is much smaller. The 
preliminary results on the same set of observed maps are con- 
sistent with SI and will be reported later. Work is in progress 
to verify them against analytical results. The full-sky maps 
and the restriction to low I < 60 (where instrumental noise is 
sub-dominant) permits the use of our analytical bias subtraction 
and error estimates. The excellent match with the results from 
numerical simulations is a strong verification of the numeri- 
cal technique. This is an important check before using Monte- 
Carlo simulations in future work targeting the / > 100 regime, 
which involves a galactic mask and a non-uniform noise matrix. 

There are strong theoretical motivations for hunting for SI 
violation in the CMB anisotropy. The possibility of non- 
trivial cosmic topology is a theoretically well-motivated pos- 
sibility that has also been observationally targeted (Ellis 1971; 
Levin 2002). Violation of SI is a generic feature of cosmic 
topology (Bond,Pogosyan & Souradeep 1998,2000). And the 
BiPS expected in flat, toroidal models of the universe has been 
computed and shown to be related to the principle directions 
in the Dirichlet domain (Hajian & Souradeep 2003a). Hence, 
the null result of BiPS have important implications for cos- 
mic topology. This approach complements the direct search 
for the signature of cosmic topology, and our results are consis- 



tent with the absence of the matched circles and the null S-map 
test of the WMAP CMB maps (Cornish et al. 1998 & 2003; 
de Oliveira-Costa et al. 1996 & 2003). Work is in progress 
to combine a full Bayesian likelihood analysis with BiPS 
analysis to constrain cosmic topology (Hajian et al. 2004), 
in particular the recently proposed dodecahedron universe 
(Luminet et al. 2003). We defer to future publication, de- 
tailed analysis and constraints on cosmic topology using null 
BiPS measurements. Other theoretical scenarios that predict 
the breakdown of SI are also being probed using BiPS, e.g., 
primordial cosmological magnetic fields (Durreretal. 1998; 
Chen et al. 2004). 

Observational artifacts such as foreground residuals, non- 
circular beam, inhomogeneous noise correlation, residual strip- 
ping patterns, etc. are potential sources of SI breakdown. Our 
null BiPS results confirm that these artifacts do not significantly 
contribute to the maps studied here. 

3. SUMMARY AND CONCLUSION 

In summary, we find null measurements of the BiPS for a se- 
lection of full-sky CMB anisotropy maps based on the first year 
of WMAP data. Our results rule out radical violation of statis- 
tical isotropy, and are consistent with null results for matched 
circles and the S-map tests of SI violation. We find that the 
excess power at low / in the best-fit theory Q with respect to 
Ci derived from the WMAP maps tends to enhance the SI vio- 
lation signatures in BiPS. Our result is an unambiguous, well- 
defined, quantitative evaluation and assessment of the SI of the 
CMB anisotropy. But pending a more careful comparison, the 
results do not necessarily conflict with a number of other sta- 
tistical tests that have suggested violation of statistical isotropy. 
We should also mention that our null BiPS results do not have 
any bearing on the non-Gaussianity of the maps. 
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Fig. 1 . — Top: Ci of the two WMAP CMB anisotropy maps. The red and magenta curves coiTespond to map A and B, respectively. The black line is a "best-fit" 
WMAP theoretical Q used for simulating SI maps. Blue dots are the average C/ recovered from 1000 realizations. Bottom: These plots show the window functions 
used. The dashed curves with increasing / coverage are "low-pass" filter, W^d,), with (, = 4, 18,40, respectively. The solid fines are "band-pass" filter W^(l,Js) 
with ft, Z,) = (13, 2), (30, 5), (30,20), (45, 20), respectively. 
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Fig. 2. — Measured values of ki for maps A and B filtered by W,''(40) (top panel) and a Gaussian filter W^^(20, 30) (bottom panel). The apparent non zero K£ at 
~ 1 cr is because the PDF is skewed to negative values and also partly because the WMAP theoretical spectrum has excess power at low / (see fig. 4). 




Fig. 3.— The probability of the two WMAP based CMB maps being SI when filtered by W/'(20,30) and a Gaussian filter W,'^(40). 
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Fig. 4. — Measured values ofBiPS for maps A and B filtered to retain power only ail = 2 and 3, assuming WMAPbf (top) and a model spectrum that matches the 
suppressed power at the lowest multipoles (Shafieloo & Souradeep 2004). 



